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Introduction 

This project focuses on the effects of weak dissipation 
on vibrational flows in migrogravity and in particular 
on (a) the generation of mean flows through viscous 
effects and their reaction on the flows themselves, and 
(b) the effects of finite group velocity and dispersion 
on the resulting dynamics in large domains. The basic 
mechanism responsible for the generation of such flows 
is nonlinear and was identified by Schlichting [21] and 
Longuct-Higgins [14]. However, only recently has it 
become possible to describe such flows self-consistcntly 
in terms of amplitude equations for the parametrically 
excited waves coupled to a mean flow equation. The 
derivation of these equations is nontrivial because the 
limit of zero viscosity is singular. This project focuses 
on various aspects of this singular problem (i.e., the limit 
C = )~ 1/2 <€ 1, where v is the kinematic viscosity 

and h is the liquid depth) in the weakly nonlinear regime. 
A number of distinct cases is identified depending on 
the values of the Bond number, the size of the nonlinear 
terms, distance above threshold and the length scales of 
interest. The theory provides a quantitative explanation 
of a number of experiments on the vibration modes of 
liquid bridges and related experiments on parametric 
excitation of capillary waves in containers of both small 
and large aspect ratio. The following is a summary of 
results obtained thus far. 

Surface-wave damping in a brimful circular 
cylinder [19] 

Henderson and Miles [6] discovered a substantial dis- 
crepancy between the measured decay rate of free os- 
cillations of a brimful circular cylinder and theoretical 
predictions based on leading order asymptotics in pow- 
ers of C 1/2 ,C = viyR* )~ x ^ 2 , where R is the radius of the 
cylinder. At the same time they found good agreement 
between the measured and predicted frequencies. The 
resolution of this discrepancy follows from the observa- 
tion by Higuera et al [7] and Martel and Knobloch [17] 
that under typical conditions the (7(C) term in the expan- 
sion of the damping rates Re( $) of capillary and gravity - 
capillary waves is comparable to the leading order term. 



Figure 1: The gravity-capillary modes for A' = 0(1). 
The dots indicate the exact solution; asymptotic results 
through 0{C ,] f 2 ) and 0(0) are indicated by dashed and 
solid lines, respectively. 

This is so despite the fact that C is typically small, 
C ^ 10 -4 . Fig. 1 shows what happens in a horizontally 
infinite layer of depth h [17]. The figure shows the de- 
cay rate Ref 5 ) and frequency Im( s) as a function of the 
wavenumber A, nondimensionalized using h 1 , for the 
parameter values used in [6] : B = pyh 2 /a = 1.96 10" 
and C = v(gh ? )~ ]/2 = 0.43 10~\ Here <7 is the surface 
tension. 

In a cylinder the corresponding calculations are 
more involved because of viscous boundary layers at 
the sidewall and the presence of a contact line. We find 

Rc(s) = C 1/2 «.'i +Cu .’2 +<7(C 3/2 ), (1) 

Im(.«) = wb-C' ,/2 wi+0(C* /a ), (2) 

where the 0{C' 12 ) terms come from viscous dissipa- 
tion in the oscillatory Stokes boundary layers near the 
solid wall and bottom of the cylinder, and the <7(C) 
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terms come from (a) viscous disipation in the bulk and 
(b) a First correction to the viscous dissipation in the 
Stokes boundary layers. The neglected 0(C 3 ^) in- 
cludes viscous dissipation in the oscillatory boundary 
layer at the free surface. We find that ^ 2/^1 for typical 
Bond numbers B = pgR 7 /<r and aspect ratios A = k/R 
is quite large even for the fundamental mode and in- 
creases with the mode number. As a result the O(C) 
term in s is important and its omission leads to errors 
that increase with the mode number, as found in [6]. 
The absence of this term in implies that the leading 
order approximation for the frequencies is much better 
than the corresponding approximation for the damping 
rates. The computation of the O(C) terms requires care 
because of a singularity in the expansion at the (pinned) 
contact line. The theoretical damping rates (1) are com- 
pared with measured rates in Table 1 , and demonstrate a 
substantial improvement over the 0(C l 12 ) results of [6]. 

The brimful cylinder is ideal for understanding the 
effects of dissipation in boundary layers because the 
meniscus is pinned to the brim, thereby eliminating un- 
known effects arising from the dynamics of the menis- 
cus. In addition the absence of comers in the container 
makes the boundary layer structure uniform and rela- 
tively simple. The damping rate calculation lends us 
confidence that we have correctly identified the major 
source of discrepancy between existing theory and ex- 
periment and that we know how to correct the theory, at 
least in cases in which surface contamination and con- 
tact angle dynamics can be ignored. It follows that (the 
bulk of the) dissipation in clean- surface experiments is 
accounted for by the classical Navier-Stokes formula- 
tion and is not due to unknown dissipation mechanisms 
in the meniscus or air-water interface. In particular any 
residual discrepancy between theory and experiment can 
now be used as a diagnostic for the presence of contami- 
nation. This fact has triggered a new experimental effort 
to measure damping rates of gravity -capillary waves in 
finite domains by M. Schatz (private communication). 
In addition the calculation forms the basis for future 
nonlinear studies of gravity-capillary waves in cylinders 
and their interactions. 

Chaotic oscillations in a nearly inviscid ax- 
isymmetric capillary bridge at 2 : 1 reso- 
nance [16] 

This paper considers a liquid bridge in microgravity sup- 
ported between two disks vibrating at two frequencies 
close to 2 : 1 resonance. Under appropriate conditions 
this vibration excites the corresponding natural vibra- 
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corresponding nondimensional damping rates, m and q 
are the azimuthal and radial wavenumbers, respectively. 
After [19]. 


lion modes (assumed axisymmetric) and these interact 
nonlinearly. In the nearly inviscid limit r < 1 this 
interaction is described by a pair of amplitude equations 
whose structure depends on the parity of the excited 
modes. Specifically, if the displacement of the disks 
follows ^ = ±A + where 

h±(i) = exp(i(fii -ffluJi )/) +.3-2 exp(t(fl 2 T ^Wi)/)], 

/< < 1, Q> = 2Qi, then the equations for the evolution 
of the amplitudes of the two competing modes take the 
form (k = 1.2) 

t&A kr = -[(1 + i)n*k\/C + n bk C - iazk((]tAk 

+ A A + — °ik$k ) T ■ * •(•!) 

Here r = fri, A = Ao C f d* -w*, t = 0( 1 ) and 

if the mode A 2 is even in N\ = A 2 Ai, N 2 = A\, This 
occurs at Ao ~ (1.249 and Ao ~ 2.23. At each of these 
values two natural oscillation modes with frequencies in 
2 : 1 ratio are present; it is these modes that are driven 
by the vibration of the two disks. Finally < measures 
the amplitude of the two modes. On the other hand, if 
A 2 is odd in ~ the nonlinearities are of third order, and 
Nk = i((a2hi | li | 2 + «-'2 1 >^2 1 2 ) -4 a- , h = 1,2. 

All the coefficients appearing in Eqs. (3) have been 
computed from the equations of motion, including the 
important (but formally subdominant) terms n bk . These 
calculations require the computation of the contribu- 
tions from the Stokes boundary layers at the disks, the 
interface boundary layer, the two comer tori near the 
edge of the disks and from the bulk. Because of the 
assumption C < 1 all nonlinear terms are purely imag- 
inary. The presence of the inhomogeneous terms in 
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Figure 2: Contour plot of the n on-oscillatory conipo- 
n cuts of the film thickness (dotted), vorticity O (solid 
for Q > 0, dashed for Q < 0) and velocity arrows for an 
ordered state. From [25]. 

Eqs. (3) implies that new phenomena are present that 
arc not captured by standard treatments of the 2 : 1 
resonance [12, 9]. Both cases have been analysed, and 
the results used to make a number of predictions about 
the sub- and superharmonic response of a liquid bridge 
subjected to this type of excitation under experimentally 
relevant conditions. The ensuing chaotic oscillations 
have beed studied numerically and characterized using 
numerically determined Liapunov exponents. 

Quasi-steady vortical structures in verti- 
cally vibrating soap films [25] 

In a recent paper Afenchenko et al [1] describe the results 
of an experimental study of flows in vertically vibrated 
soap films; in this experiment the frame holding the film 
was rigidly and symmetrically attached to the lateral wall 
of a closed cavity, which was vibrated vertically. These 
vortical structures (first reported by Taylor [24]) are 
similar to ones described by Airiau [2], who excited 
the film using a loudspeaker fitted to the bottom of the 
cavity, which was otherwise open. We have examined 
in detail the possible mechanisms that could lead to 
the generation of such flows, with special emphasis on 
those effects that are independent of the shape of the 
frame supporting the film, of the attachment mode (i.e., 
the size of the meniscus), and of the excitation process. 


X 

Figure 3*. Vorticity contour* for a spatially disordered 
state. From [25]. 

Thus we focus on those effects that arc present in any 
experiment, e.g., in [1, 2]. In this theory the observed 
vortical structures are a consequence of the oscillatory 
tangential and normal stresses on the film due to the 
surrounding air. The air also damps out Marangoni 
waves. Because of nonlinearities these stresses (and the 
inertia of the film) produce non-oscillatory deflection 
of the the film, variation in its thickness, and forcing 
terms that combine to produce a streaming flow in the 
film. Coupled evolution equations for these quantities 
are derived. When the excitation frequency is close 
to an eigenfrcquency of a Marangoni mode of the soap 
film, the dominant forcing in these equations arises from 
the nonlinear hydrodynamics within the film volume, 
while both volume forcing and surface forcing by air arc 
important when there is no resonance with a Marangoni 
mode. The computed vortex patterns (see Figs. 2,3) 
agree qualitatively with the experiments. 

Compressional modes in parametrically driven 
Faraday waves in an annulus [18] 

Experiments by Douady et al [4] on parametrically 
driven water waves in an annulus (C = 4.4 10~ 4 , 

B = 8.9) reveal the presence of a secondary insta- 
bility of a uniform pattern of standing waves (SW) in 
the form of an oscillatory compression mode (CM). We 
develop a theory describing this observation based on 
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(suitably modified) nonlocal amplitude equations for the 
subharmonic response of the system of the type derived 
by Knobloch and De Luca [10] for traveling wave con- 
vection and Alvarez-Pereira and Vega [3] for pulsating 
flames, with the nonlinear coefficients deduced from 
Hansen and Alstrom [5]: 

,4 1 = (-1 + iv)A + fi{B) + + 7<|B| 2 )) + ioA, x 

(4) 

B, = (— 1 + 1 v)B + /f (A) + iB(3\B\“ 4- 7'{M| 2 )) +i<*B xx , 

(5) 

subject to the periodic boundary conditions A{x -f 1 J) = 
B(x + 1, /) = Here (* • •) represents a 

spatial average over 0 < x < 1. These equations de- 
scribe the slow temporal evolution of the amplitudes of 
lcft-(. l) and right-traveling (B) waves in their comoving 
frames, and provide an asymptotic description of the 
system sufficiently close to threshold of the Faraday 
instability that the envelope dynamics are dominated 
by advection at the group velocity of the waves. A 
description of this type is appropriate in the weak dissi- 
pation limit in which the natural modes of the unforced 
system, left and right-traveling waves, decay slowly. 
For this reason the coefficients of the cubic terms are 
purely imaginary. In the first instance we ignore the 
presence of viscosity-generated mean flows and study 
the resulting nonlocal equations. Within this description 
the onset of the compression mode is described as a sec- 
ondary pitchfork bifurcation from SW with eigenvector 
that breaks the SW reflection symmetry. We show 
that such CM are described by steady but nonuniform 
solutions of the nonlocal equations with periodic bound- 
ary conditions (appropriate to an annular cell) and are 
present for experimental parameter values. We explore 
the transition from this mode and identify consecutive 
bifurcations to nonsteady but uniform states, and to 
states with complex spatio-temporal dynamics (Fig. 4), 
some of which are strongly hysleretic. Future work will 
include viscosity-driven mean flows with self-consistent 
coupling to the amplitude equations and explore the dy- 
namics of the resulting system of equations. Such mean 
flows arc present not only in large aspect ratio annuli 
but also in smaller systems, and hence are relevant to 
the experiments. 

A related investigation of the corresponding prob- 
lem with exact or broken Da symmetry is under way, 
motivated by the Faraday problem in a square or nearly 
square container. In a container with exact square sym- 
metry two modes related by a 90° rotation are exactly 
degenerate; this degeneracy is broken if the container is 
slightly rectangular. Simonelli and Gollub [22] found 



Figure 4: Stable solutions of Eqs. (4,5) for v = 0, 3 = 3, 
y = —1, a = 0.1, showing SW (0 < // < 1.46), CM 
(1.46 < (,i < 1.62) and complex states (// > 1.62). 


no time-dependent dynamics associated with the result- 
ing mode interaction in a square domain. However in 
a slightly rectangular container they observed periodic 
and chaotic bursts very close to onset. The real parts 
of the necessary cubic coefficients are being calculated 
to leading order in 0{C) together with the mean flows 
generated in these two geometries. A self-consistent 
inclusion of such flows offers the possibility of provid- 
ing the first quantitative description of this interesting 
behavior. An explanation along the lines described in 
the next section for the presence of bursts is anticipated. 

Bursts [15] 

Under appropriate conditions the competition between 
nearly degenerate Hopf modes with odd and even parity 
results in dramatic bursting behavior very near the onset 
of primary instability. This behavior arises in two- 
dimensional systems of large but finite aspect ratio 
undergoing an oscillatory instability, or in systems with 
nearly square symmetry, as discussed by Landsberg and 
Knobloch [13]. In both cases it is due to the same 
mechanism. 

In a slender system with left-right reflection sym- 
metry (such as a narrow rectangular convection cell) 
undergoing an oscillatory instability from the trivial 
state the first two unstable modes typically have op- 
posite parity under reflection. Moreover, because the 
neutral stability curve for the unbounded system has 
a parabolic minimum these set in in close succession 
as the control parameter is increased. We write the 
perturbation from the trivial state as 

y,i) = r.* Re {z+f+{x. y) + z-f-(x,y)} + <f>(>), 

where f. < 1, = ±f±(x,y), and y denotes 

the transverse variables. The complex amplitudes ~±(*) 
then satisfy the equations [ 1 3] 

i± = [A± AA + i(-’± A^)]s±+.4(|; + | 2 + |^-| 2 );± 

+ 5 |;±| 2 ;±+<^± 4 - ( 6 ) 
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In these equations the nonlinear terms have identical 
(complex) coefficients because of an approximate inter- 
change symmetry between the odd and even modes. The 
resulting Da symmetry [13] is weakly broken whenever 
A A / 0 and/or A^ ^ 0, a consequence of the finite 
aspect ratio of the system. The dynamics in systems 
with exact Da symmetry (A A = Au; — 0) can be shown 
to take place on a two-dimensional manifold [23] and 
hence are necessarily simple although the solutions can 
become unbounded. This is no longer so when the 
Da symmetry is broken. In this case the solutions can 
experience repeated episodes of dramatic growth fol- 
lowed by collapse but do not become unbounded. The 
resulting bursts can be periodic or chaotic and set in at 
a secondary instability very close to onset. To identify 
such bursts we write 


2 + 7 ± 4 


f . i (± 0 + t ')/2 


and introduce a new time-like variable r defined by 
dr /dt = p ~ l . In terms of these variables Eqs. (6) 
become 


dp 

dr 


-p[i\ i{ + Br( 1 + cos 2 0) + C R sin 2 0 cos 2<p] 




— 2( A 4- AA co±Q)p 


(I) 


Figure 5: Bursts arising from periodic and chaotic ro- 
tations. (a) A = 0.1 and (b) A = 0.072. The coefficients 
are AA = 0.03, A^ = 0.02, A = 1 -1.5/, B = -2.8 + 5f, 
C = 1 + i. 


- — = sin 0 [cos 0( — Br -f Cr cos 2<;6) — Cj sin 2<p] 

dr 

— 2AA p sin 0 (5) 

— — = cos 0(B i — Ci cos 2q) — Cr sin 2(f> 

dr 

+2Au; p. (9) 

where. I = Ar + iAj, etc., together with a decoupled 
equation for f (/). The amplitude of the disturbance 
is measured by r = \z+\ 2 H- \z-\ 2 = p” 1 ; thus p = 0 
corresponds to infinite amplitude states. Eqs (7-9) show 
that the restriction to the invariant subspacc E = {p = 0} 
is equivalent to taking AA = Avt? = 0 in (8,9). The 
fixed points of the resulting D 4 -symmetric problem 
correspond to (infinite amplitude) periodic oscillations 
in time because of the decoupled phase f (f). Depending 
on A, B and C the subspace E may contain additional 
fixed points and/or limit cycles [23]. In our scenario, 
a burst occurs for A > 0 when a trajectory follows 
the stable manifold of a fixed point (or a limit cycle) 
Pi e E that is unstable within E. The instability within 
E then kicks the trajectory towards another fixed point 
(or limit cycle) P 2 € E. If this point has an unstable 
p eigenvalue the trajectory escapes from E towards a 
p > 0 fixed point (or limit cycle), forming a burst. 
If AA and/or Au; ^ 0 this fixed point may itself be 
unstable to perturbations of type P\ and the process then 


repeats. The scenario thus requires that at least one of 
the branches in the /^-symmetric system be subcritical 
(Pi) and one supercritical (P 2 ). 

When AA and/or Aa.< ^ 0 two types of oscillations 
in are possible: rotations and librations. These 
oscillations are coupled to excursions in amplitude. 
Fig. 5 shows a typical sequence of large amplitude bursts 
arising from repeated excursions towards the infinite 
amplitude (p = 0) solutions. The amplitude of these 
events decreases with increasing A A and their frequency 
increases, much as found in several experiments. The 
mechanism outlined here arises naturally in systems 
with reflection symmetry and appears to be responsible 
for the regular and irregular bursting exhibited by such 
systems very close to onset of a primary oscillatory 
instability. 

Dynamics of parametrically modulated dis- 
sipative systems in an annulus [11] 

Dynamics of parametrically modulated dissipative sys- 
tems undergoing a symmetry -breaking Hop f bifurcation 
on a line are explored in full generality, with particular 
emphasis on the case in which either standing waves or 
travelling waves are subcritical [20]. The spatially uni- 
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form states satisfy equations of the form (cf. Eqs. (4,5)) 
A = a A + + c\A\ 2 A + d]B\ 2 A (10) 

D = aB + bA + c\B\ 7 B + d\A\ 2 B. ( 11 ) 

Here a, <■ and d are complex while h may be taken to 
be positive. In contrast to Eqs. (4,5) we assume that 
Re(«) can pass through zero leading to a spontaneous 
oscillatory instability, and study the effect of parametric 
modulation (b > 0) on these oscillations. In closely 
related equations (see Eqs. (6)) bursts occur when one 
of the branches is subcritical (cf. [8]). We have identi- 
fied the possible dynamical regimes involving spatially 
uniform slates. These include standing waves and non- 
symmetric mixed modes phase-locked to the drive and 
interactions between them. The effects of group velocity 
and dispersion on the stability properties of the spatially 
uniform states identified above are being analyzed. As 
a result a number of new instabilities leading, for ex- 
ample, to compression-like states have been described 
and these lead to various novel spatially inhomogeneous 
stales. 

Acknowledgements 

This work was supported by NASA under grant UGS97- 
0308 and by NSF under grant DMS-9703684. We are 
grateful to F. Higuera, M. Higuera, F. Manccbo, C. 
Martel, J. Moehlis and J. Nicolas for their assistance 
in developing the projects described above, and to P. 
Weidman for making his experimental results available 
prior to publication. 

References 

[1] V.O. Afenchenko, A.B. Ezersky, S.V. Kiyashko, 
M.I. Rabinovich and P.D. Weidman, Phys. Fluids, 
in press. 

[2] M. Airiau, DEA Report, Ecole Normale Su- 
pcricure, Paris (1986). 

[3] C. Alvarez- Pereira and J.M. Vega, Euro. Jnl. of 
Appl. Math. 3, 55 (1992). 

[4] S. Douady, S. Fauve and O. Thual, Europhys. Lett. 
10,309 (1989). 

[5] P.L. Hansen and P. Alstrom, J. Fluid Mech. 351 , 
301 (1997). 

[6] D.M. Henderson and J.W. Miles, J. Fluid Mech. 
275 , 285 (1994). 


[7] M. Higuera, J.A. Nicolas and J.M. Vega, Phys. 
Fluids 6, 438 (1994). 

[8] P.C. Hirschberg and E. Knobloch, Physica D 90, 
56(1996). 

[9] D.W. Hughes and M.R.E. Proctor, J. Fluid Mech. 
244, 583(1992). 

[10] E. Knobloch and J. DeLuca, Nonlinearity 3, 575 
(1990). 

[11] E. Knobloch, C. Martel, J. Moehlis and J.M. Vega, 
in preparation. 

[12] E. Knobloch and M.R.E. Proctor, Proc. R. Soc. 
London A 415 , 61 (1988). 

[13] A.S. Landsberg andE. Knobloch, Phys. Rev. E 53, 
3579(1996). 

[14] M.S. Longuet-Higgins, Phil. Trans. R. Soc. Lon- 
don A 245, 535 (1953). 

[15] J. Moehlis and E. Knobloch, Phys. Rev. Lett., in 
press. 

[16] F.J. Mancebo, J.A. Nicolas and J.M. Vega, Phys. 
Fluids, in press. 

[17] C. Martel and E. Knobloch, Phys. Rev. E 56, 5544 
(1997). 

[18] C. Martel, E. Knobloch and J.M. Vega, in prepa- 
ration. 

[19] C. Martel, J.A. Nicolas and J.M. Vega, J. Fluid 
Mech. 360,213(1998). 

[20] H. Riecke, J.D. Crawford and E. Knobloch, in 
NATO ASI Series B 237, 61 (1991). 

[21] H. Schlichting, Z. Phys. 33, 327 (1932). 

[22] F. Simonelli and J.P. Gollub, J. Fluid Mech. 199 , 
471 (1989). 

[23] J.W. Swift, Nonlinearity i, 333 (1988). 

[24] S. Taylor, Proc. R. Soc. London 27, 71 (1878). 

[25] J.M. Vega, F.J. Higuera and P.D. Weidman, J. 
Fluid Mech., in press. 


482 


